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ABSTRACT 

We measure the 3D genus topology of large scale structure using Luminous Red Galaxies (LRGs) in 
the Sloan Digital Sky Survey and find it consistent with the Gaussian random phase initial conditions 
expected from the simplest scenarios of inflation. This studies 3D topology on the largest scales ever 
obtained. The topology is sponge-like. We measure topology in two volume-limited samples: a dense 
shallow sample studied with smoothing length of 21/i _1 Mpc, and a sparse deep sample studied with a 
smoothing length of 34/i _1 Mpc. The amplitude of the genus curve is measured with 4% uncertainty. 
Small distortions in the genus curve expected from non-linear biasing and gravitational effects are well 
explained (to about 1-a accuracy) by ./V-body simulations using a subhalo-finding technique to locate 
LRGs. This suggests the formation of LRGs is a clean problem that can be modeled well without any 
free fitting parameters. This bodes well for using LRGs to measure the characteristic scales such as 
the baryon oscillation scale in future deep redshift surveys. 

Subject headings: large-scale structure of universe - cosmology: observations - methods: numerical 
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1. INTRODUCTION 

Topology of large-scale structure in the universe has 
been studied over the years through analyses of large- 
scale distribution of galaxies in three dimensions. The 
genus statistic was used to characterize the topology 
(Gott, Melott, & Dickinson 1986; Hamilton, Gott, & 
Weinberg 1986; Gott, Weinberg, & Melott 1987; Gott, et. 
al. 1989; Vogeley, et. al. 1994; Park, Kim, & Gott 2005a; 
Park et. al. 20055). Measuring the genus as a function 
of density (Gott et al. 1986, 1987) allows us to com- 
pare the topology observed with that expected for Gaus- 
sian random phase initial conditions, as predicted in a 
standard big bang inflationary model (Guth 1981; Linde 
1983) where structure originates from random quantum 
fluctuations in the early universe. 

At large scales where fluctuations are still in the lin- 
ear regime, growing in place without changing topology, 
topology at the median density level is expected to be 
sponge-like (Gott et al. 1986, 1987). Small deviations 
from the random phase curve give important informa- 
tion about biased galaxy formation and non-linear grav- 
itational clustering as shown by perturbation theories 
and large iV-body simulations (Matsubara 1994; Park, 
et al. 2005a). Previous models of galaxy clustering 
suggested either a meatball topology, isolated clusters 
growing in a low density connected background (Press 
& Schechter 1974; Soneira & Peebles 1978), or a Swiss 
cheese topology, isolated voids surrounded on all sides 
by walls (Einasto, Joeveer, & Saar 1980). But studies of 
many observational samples have shown in every case a 
sponge-like median density contour as expected from in- 
flation (Gott et al. 1986, 1989, 2008; Moore et al. 1992; 
Vogeley et al. 1994; Canaveses et al. 1998; Hikage et al. 
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2002, 2003; Park et al. 20056). 

We have previously measured the 3D topology of the 
nearby Main galaxies (Strauss et al. 2001) of the Sloan 
Digital Sky Survey (York et al. 2000; Stoughton et al. 
2002; Gunn et al. 2006) data at smoothing lengths rang- 
ing from 3.5 to 11 h Mpc (see Park et al. 20056 for 
details). The overall agreement with the best-fit Gaus- 
sian random phase genus curve and the mock surveys in 
a ACDM is remarkable, strongly supporting the predic- 
tions of inflation. Small deviations from the Gaussian 
random phase distribution are expected because of non- 
linear gravitational evolution and biased galaxy forma- 
tion, and these can now be observed with sufficient accu- 
racy to do model testing of galaxy formation scenarios. 
In this paper we will analyze the 3D distribution of the 
luminous red galaxys (LRGs; Eisenstein et al. 2001) ob- 
served by the SDSS to study the genus topology at very 
large scales. The results will be compared with those 
from the nearby galaxies and with those from N-body 
simulations of the inflation-based ACDM universe. 

2. LUMINOUS RED GALAXIES 

The LRG galaxies in the SDSS have already proven 
useful in refining cosmological parameters (Tegmark et 
al. 2006). Measuring the genus topology of the LRGs is 
of particular importance because it allows us to observe 
topology on the largest scales, and because the LRGs are 
expected to play an important role in characterizing dark 
energy, i.e. through the ratio of dark energy pressure to 
energy density (cf. Bessett et al. 2005; Eisenstein et al. 
2005; and Perceival et al. 2007 for a review). 

We use the LRG sample from the SDSS. The sky cov- 
erage of the sample is similar to that of SDSS data Re- 
lease 5 (Adelman-McCarthy et al. 2006). The angular 
selection function is based on the spherical polygon de- 
scription of LSS-DR4plus catalog (Blanton et al. 2003, 
2005). To maximize the volume-to-surface ratio we trim 
the DR4plus sample (see Figure 1 and the text of Park 
et al. [2007]). These cuts leave 68,203 LRGs over about 
4464 deg 2 in the survey region. 

We constructed two volume-limited subsamples of 
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LRGs in passively evolved luminosity and redshift 
ranges. The subsamples are —21.2 < M ffl o.3 < —23.2 
with 0.2 < z < 0.36, and -21.8 < M s>0 .3 < -23.2 with 
0.2 < z < 0.44, where we have applied fc-corrections 
and passively evolved the galaxies to a fiducial redshift 
of 0.3. We label them "SHALLOW" and "DEEP", re- 
spectively. Both subsamples are selected to reduce the 
artifacts caused by the LRG selection cuts. The result- 
ing comoving number density of LRGs are approximately 
constant, (see Zehavi et al. 2005 for details). The sam- 
ple DEEP (SHALLOW) contains 32,067 (14,968) galax- 
ies. The absolute magnitude and redshift limits of these 
samples are illustrated in Figure 1. We will measure the 
genus topology of the DEEP (SHALLOW) sample with 
a smoothing length of 21 (34) ft, _1 Mpc. The smoothing 
lengths chosen are large enough to avoid unwanted shot 
noise effects (c.f. Gott, Melott & Dickinson 1985) 
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3. N-BODY SIMULATION 

We completed a 2048 3 particle 4915.2/i~ 1 Mpc size Cold 
Dark Matter Simulation adopting the WMAP 3-year 
parameters (Spergel et al. 2007); f2 m = 0.238, Q\ = 
0.762, n b = 0.042, n s = 0.958, a s = 0.761, and h = 0.732, 
where fi m , &b are the density parameters due to mat- 
ter, cosmological constant, and baryon, respectively. n s 
is the slope of the power spectrum at large scales, and 
og is the rms density fluctuation at 8/i _1 Mpc scale. The 
simulation is run by a PM+Tree N-body code (Dubinski 
et al. 2004) and has a force resolution of 0.24 h Mpc. 

We identify the most massive dark matter subhalos as 
LRGs. The subhalos are the gravitationally bound ha- 
los not subject to tidal disruption (Kim & Park 2006), 
and having more than 30 member particles. To com- 
pare with the observational data we saved the parti- 
cle positions and velocities along the past light cone 
of 12 well-separated observers who are making mock 
SDSS DR4plus LRG surveys. The subhalos are found on 
the past light cone surface. The minimum halo masses 
matching the observed number densities of the SHAL- 
LOW and DEEP LRG samples are found as a function of 
redshift and applied to the subhalos data on the past line 
cone when mock surveys are made. In each mock survey 
we pick the most massive subhalos when the mean num- 
ber density of the LRGs is matched. At small non-linear 
scales the genus depends sensitively on the method by 
which galaxies are identified within the simulation. Park 
et al. (2005a) and Gott et al. (2008) have extensive 
studies on this question. 

4. THE GENUS AND RELATED STATISTICS 

To measure the genus we smooth the LRG number 
density distribution with a Gaussian smoothing ball of 
radius Rg- To remove the small radial variation of LRG 
number density in our samples, we divide the density by 
the radial selection function in the Appendix of Zehavi 
et al. (2005). Then the iso-density contour surfaces of 
the smoothed galaxy density distribution are searched to 
calculate the genus. The contour surfaces are labeled by 
is, which is related with the volume fraction / on the high 
density side of the density contour surface by 
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Fig. 1. — Distribution of Luminous Red Galaxies of the SDSS. 
Solid boundary lines define the SHALLOW sample, and the dashed 
lines define the DEEP sample. 

The / = 50% contour corresponds to the median volume 
fraction contour (y = 0). For Gaussian random phase 
initial conditions the genus curve is: 



7(^=A(l-^)e- 2 / 2 , 



(2) 



where the amplitude A = ((fc 2 )/3) 3/2 /27r 2 and (fc 2 ) is 
the average value of k 2 in the smooth power spectrum 
(Hamilton, Gott, & Weinberg 1986; Doroshkevich 1970) 
The shape of the genus curve can be parameterized 
by several variables. First, there is the amplitude of the 
genus curve as measured by fitting the amplitude of the 
best fitting Gaussian random phase curve from Equation 
2. This gives information about the power spectrum and 
phase correlation of the density fluctuation. Deviations 
of the genus curve from the theoretical random phase 
case can be quantified by the following three variables: 



Av 



\\g(v)vdv 
J-ifl'rfM du ' 



A v = 



1-2.2 ffrf M d " 



Ac 



/l.2 ffrfM dp ' 



(3) 



(4) 



(1) 



where g?i{y) is the genus of the best-fit random phase 
curve (Eq. 2). Thus Av measures any shift in the cen- 
tral portion of the genus curve. The Gaussian curve (Eq. 
2) has Av = 0. A negative value of Av is called a "meat- 
ball shift" caused by a greater prominence of isolated 
connected high-density structures which push the genus 
curve to the left. Ay (Ac) measure the observed num- 
ber of voids (clusters) relative to those expected from the 
best-fitting Gaussian curve. 

5. RESULTS 

We have measured the 3D genus topology of the LRGs 
in the SHALLOW and DEEP samples with the smooth- 
ing lengths described above. The genus curves are shown 
in Figure 2. The maximum value of the genus curve 
reaches about 170, and the amplitude is measured with 
only 4% error at 21/i _1 Mpc scale. This compares with 
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Fig. 2. — Genus curves for observed (two noisy curves) and sim- 
ulated (two long-dashed curves) Luminous Red Galaxies. Gaus- 
sian fits (solid lines) are also shown. The genus curves with 
higher amplitudes are for the SHALLOW sample smoothed with 
Rg = 21h~ 1 TsApc, and those with lower amplitudes are for the 
DEEP sample with Rq = 34/t- 1 Mpc. 



the previously smallest error of 9.4% at 5 ft^Mpc scale 
for the genus obtained from the Main Galaxy Sample 
of the SDSS (Park et al. 20056). The data for each 
sample is shown as the jagged line. Each genus curve is 
compared with the best-fit Gaussian random phase curve 
from Equation 2 (smooth solid lines). The observational 
data agree well with the Gaussian curves, but there are 
small differences between them due to the non-linear bi- 
asing and gravitational effects (Park et al. 2005a). 

The dashed lines show the mean of 12 mock SDSS Sur- 
veys constructed from our 2048 3 particle simulation. Al- 
though our two volume-limited samples contain galax- 
ies typically with different absolute magnitude and red- 
shift, the genus curves measured from them show a very 
good overall agreement. Moreover, there is a remark- 
able agreement between observations and the mock sur- 
veys. The number of voids (the depth of the valley at 
v < —1), is less than the number of clusters (the depth 
of the valley at v > 1). The phenomenon that clusters 
outnumber voids is reproduced by the simulations ex- 
tremely well. The amplitude of the genus curve agrees 
well between observations and simulations, showing that 
the non- linear power spectrum and degree of phase corre- 
lation of the LRG distribution smoothed at Rq = 21ft. -1 
Mpc and 34/i _1 Mpc are also modeled well. The small 
deviations of the genus curves from the Gaussian cases 
are quantified by our distortion parameters Ais, Ay, and 
Ac- These statistics for the observed LRG distribution 
are within approximately la of those expected from the 
12 mock catalogs as shown in Figures 3 and 4. 

In Figure 3 the predictions for random phase curve 
(Eq. 2) are shown by stars with three and five legs for 
the cosmological parameters estimated from WMAP 3- 
year (Spergel et al. 2007) and 5-year maps only (Dunklcy 
et al. 2008), respectively. In Figure 4 the random phase 
curve is indicated by a point where Ay = 1 and Ac = 1. 
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Fig. 3. — Genus per unit volume multiplied by the cube of the 
smoothing length: gRc 3 versus the shift parameter Ais. Filled cir- 
cles with error bars are for the SDSS LRG galaxies in the SHAL- 
LOW sample (marked as LRG z m ax = 0.36) and in the DEEP 
sample. Open circles are the means of 12 mock SHALLOW and 
DEEP LRG surveys performed along the past light cone surface in 
a ACDM universe. The filled square is from the SDSS galaxy dis- 
tribution at 6 /i -1 Mpc scale taken from Park et al. (2005a). The 
stars with 3 and 5 legs are the linear theory predictions for the 
WMAP 3 year and 5 years cosmological parameters, respectively. 
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Fig. 4. — Same as Fig. 3 but for the Ac and Ay parameters. The 
cross on the right is the prediction for a random phase fluctuation. 



The mock catalogs are indicated by open circles. The 
observational data are indicated by filled circles with la 
error bars showing the la variation from mock survey to 
mock survey. Importantly, there is no free fitting param- 
eter here for any of these quantities. So it is noteworthy 
that in A, Ais, Ay, and Ac the data are within about la 
of the predicted values. The statistic showing the max- 
imum difference between observation and simulation is 
Ay; the observed voids are somewhat larger and fewer 
relative to those found in mock SDSS LRG surveys. 

This suggests that finding the objects corresponding to 
LRGs is a relatively clean problem. Our genus topology 
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TABLE 1 
Genus-Related Statistics 



Name 


Magnitude 




Redshift 


d 


Rg 


s 


Amplitude 


Av 


A v 


Ac 


SHALLOW 
DEEP 


-21.2 < M g ,o.3 < 
-21.8 < Afg.o.3 < 


to to 

CO CO 

to to 


0.2 < z < 0.36 
0.2 < z < 0.44 


22.1 
35.0 


21 
34 


7 
11 


167.4 ± 7.0 
79.6 ±6.0 


-0.041 ±0.026 
-0.056 ± 0.030 


0.694 ± 0.065 
0.801 ±0.070 


1.087 ±0.057 
1.127 ±0.092 



Note. — Cols. (1) Name of samples, (2) Rest-frame g-band absolute magnitudes, (3) Redshift, (4) Mean separation of galaxies in units 
of h~ 1 Mpc, (5) Smoothing length in units of h~ 1 Mpc, (6) pixel scale in units of h~ 1 M.pc, (7) The amplitude of the observed genus curve, 
(8) The shift parameter, (9) and (10) are the cluster and void abundance parameters, respectively. Uncertainty limits are estimated from 
mock surveys in redshift space. 



study shows that iV-body simulation together with iden- 
tification of LRGs with the most massive dark subhalos 
can produce robust results in very good agreement with 
the observational data. Our results here show that the 
ACDM model we adopted in the iV-body simulation is a 
close approximation of the real universe. The results also 
suggest that it will be possible to model well the forma- 
tion of LRGs in the future Sloan III survey by using large 
A^-body simulations to check and correct for any small 
systematic effects due to non-linear gravitational effects 
and biasing without needing free fitting parameters. 

6. CONCLUSIONS 

We measured the 3D genus topology of large scale 
structure using the LRG galaxies in the Sloan Survey. 
The amplitude of the genus curve is now measured within 
about 4% error at 21 /i -1 Mpc smoothing scale. The re- 
sults are consistent with the predictions of the ACDM 
model with WMAP cosmological parameters, and with 
Gaussian random phase initial conditions expected from 
standard inflation. The topology is sponge-like. Small 
distortions from the random phase curve expected from 
non-linear gravitational effects and biasing are observed. 
These are modeled well (at the la level) by A^-body sim- 
ulations that pick the largest gravitationally bound halos 
not subject to tidal disruption as the sites for LRG galax- 
ies to form. The data on topology, including the number 



of voids and clusters, the location of the peak in the genus 
curve and its amplitude are all fit well without any free 
fitting parameters. Thus, the formation of LRGs appears 
to be a relatively clean problem that can be modeled well 
using cold dark matter and gravity alone, and this bodes 
well for using A^-body simulations to calibrate measure- 
ment of baryon oscillations using LRGs. 
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